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ABSTRACT 

Diffuse interstellar bands (DIBs) have been discovered for almost a century, but their 
nature remains one of the most challenging problems in astronomical spectroscopy. 
^ ^ Most recent work to identify and investigate the properties and carriers of DIBs con- 

centrates on high-resolution spectroscopy of selected sight-lines. In this paper, we 
report detections of DIBs in the Sloan Digital Sky Survey (SDSS) low-resolution spec- 
tra of a large sample of Galactic stars. Using a template subtraction method, we have 
successfully identified the DIBs AA5780, 6283 in the SDSS spectra of a sample of about 
2,000 stars and measured their strengths and radial velocities. The sample is by far 
the largest ever assembled. The targets span a large range of reddening, Eb-v ^ 0.2 
- 1.0, and are distributed over a large sky area and involve a wide range of stellar pa- 
rameters ( effective temperature, surface gravity and metallicity), confirming that the 
carriers of DIBs are ubiquitous in the diffuse interstellar medium (ISM) . The sample 
is used to investigate relations between strengths of DIBs and magnitudes of line-of- 
sight extinction, yielding results (i.e., EW{5780) = 0.61 x Eb-v and £;W^(6283) = 
1.26 X Eb-v) consistent with previous studies. DIB features have also been detected 
in the commissioning spectra of the Guoshoujing Telescope (LAMOST) of resolving 
power similar to that of SDSS. Detections of DIBs towards hundreds of thousands of 
stars are expected from the on-going and up-coming large scale spectroscopic surveys 
such as RAVE, SDSS III and LAMOST, particularly from the LAMOST Digital Sky 
Survey of the Galactic Anti-center (DSS-GAC). Such a huge database will provide 
an unprecedented opportunity to study the demographical distribution and nature of 
DIBs as well as using DIBs to probe the distribution and properties of the ISM and 
the dust extinction. 

Key words: ISM: lines and bands - ISM: dust, extinction - stars: general. 



1 INTRODUCTION 

Diffuse interstellar bands (DIBs) are absorption features 
ubiquitously detected in the spectra of reddened stars. They 
are observed from the near UV, optical, to the near infrared, 
and arise from the interstellar medium (ISM). DIBs were 
first discovered by Heger (1922). To date over 400 DIBs 
have been detected (Sarre 2006; Hobbs et al. 2008, 2009). 
DIBs are also detected in extragalactic sources, such as the 
Small Magellanic Cloud (Ehrenfreund et al. 2002; Welty 
et al. 2006; Cox et al. 2007), the Large Magellanic Cloud 
(Ehrenfreund et al. 2002; Welty et al. 2006; Cox et al. 2006), 
the spiral galaxy NGC 1448 (SoUerman et al. 2005), the An- 



dromeda galaxy M31 (Cordiner et al. 2008a, 2011), the Tri- 
angulum galaxy M33 (Cordiner et al. 2008b), dusty star- 
burst galaxies (Heckman & Lehnert 2000) and Ca II and 
damped Lya absorbers towards quasars at cosmological dis- 
tances (Junkkarinen et al. 2004; York et al. 2006; Ellison et 
al. 2008; Lawton et al. 2008). 

The most prominent DIBs in the optical, e.g., AA4429, 
5780, 57970, are well studied. They are known empirically to 
trace the neutral hydrogen phase of the ISM, with intensities 
that correlate well with the line-of-sight color excess Eb-v, 
the neutral hydrogen column density and the Na I column 
density (e.g., Herbig 1993; Friedman et al. 2011). Correla- 
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^ We cite all DIBs with reference to their nominal air wave- 
lengths. 
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tions between different DIBs deduced from samples of spec- 
tra towards various sight-lines have been used to infer the 
origin(s) of the bands - sets of well-correlated DIBs probably 
arise from the same or similar carriers (Krelowsky & Walker 
1987; Josafatsson &: Snow 1987; Westerlund & Krelowsky 
1989; Cami et al. 1997; Moutou et al. 1999; Wszolek & God- 
lowski 2003; Friedman et al. 2011). Although DIBs have been 
widely detected and studied, none of them have been spec- 
troscopically identified and their carriers (widely believed to 
be grains or molecules) remain elusive (Herbig 1995; Sarre 
2006). 

DIBs exhibit a wide range of widths, with full widths at 

half maximum (FWHM) ranging from below 1 A to tons of 
A. High-resolution spectroscopy of reddened OB stars has 
been used to search for new features and to study the intrin- 
sic profiles, fine structures and possible relations with other 
sharp interstellar absorption lines (e.g., Sarre et al. 1995; 
Galazutdinov et al. 2003; Thorburn et al. 2003; Hobbs et al. 
2008, 2009; McCall et al. 2010; Friedman et al. 2011; Vos 
et al. 2011). Luminous OB stars arc targeted because their 
spectra possess few photospheric absorption lines, making it 
much easier to detect and measure the relatively weak DIBs. 
Hitherto, more than one hundred high-resolution spectra of 
OB stars have been accumulated. For DIBs of intermedi- 
ate widths, such as the strong DIBs AA5780, 6283, low-to- 
intermediate resolution spectroscopy matches their widths 
very well, and thus can be used to detect and study them 
efficiently. At a resolving power of ~ 3,000, Cordiner et al. 
(2008a, 2008b, 2011) detected the DIBs AA5780, 6283 and 
several other features in the optical spectra of a number 
of supergiants in M33 and M31 with the DEIMOS spec- 
trograph of the Keck telescope. Using the Radial Veloc- 
ity Experiment (RAVE, Steinmetz et al., 2006; Zwitter et 
al. 2008; Siebert et al. 2011) data of a resolving power of 
~ 7,500, Munari et al. (2008) detected the DIB A8621 to- 
wards 68 hot stars. They found a tight correlation between 
strengths of the DIB A8621 and reddening, implying that 
the feature could bo a good measure of extinction. This fear 
ture will be targeted by the forth-coming next generation 
astrometric satellite GAIA at a resolving power of ~ 11,500 
(Katz et al. 2004). Using blue- violet spectra at a resolving 
power of 2,500 from the Galactic 0-star spectroscopic survey 
(GOSSS, Mai'z Apellaniz et al. 2011), Penades Ordaz et al. 
(2011) detected DIBs AA4501,4726,4762,4780,4964 towards 
a few hundred O stars. 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) 
has obtained low-resolution (A A 3.3 A at 6,000 A) spec- 
tra of over 500,000 stars in its Data Release 8 (DR8; Aihara 
et al. 2011). A fraction of the stars are towards high extinc- 
tion sight-lines. In this paper, we report detections of DIBs 
in the SDSS spectra of a variety of different types of stars. 
DIBs AA5780, 6283, two strong DIB features of intermedi- 
ate widths that are well matched with the SDSS spectral 
resolution, are selected for the investigation. 

The paper is organized as follows. In Section 2, we in- 
troduce the data and method used to identify and mea- 
sure DIBs. The results of DIBs AA5780, 6283 detected in 
the SDSS spectra are presented in Section 3. Detections of 
DIB features in the commissioning spectra of the Guosliou- 
jing Telescope (LAMOST, Wang et al. 1996) are reported in 
Section 4. The prospect and implications of detecting and 
studying DIBs towards a huge number of sight-lines from 



several on-going and up-coming large scale low-resolution 
spectroscopic surveys are discussed in Section 5. 



2 DATA AND METHOD 

We use a template subtraction method to detect DIBs in 
the SDSS spectra. This method recjuires a sample of (red- 
dened) target stars and an (unreddened/low- reddening) con- 
trol sample of stars of similar spectral type. Target stars are 
defined as stars whose spectra are expected to display DIB 
features. Control stars are those whose spectra posses ex- 
tremely weak or no DIB features, so that the spectra can 
be used to construct spectral templates containing no DIBs. 
In this Section, selections of target and control stars from 
the SDSS spectroscopic survey as well as the technique of 
template subtraction are introduced. 



2.1 Data 

Both samples of target and control stars are selected from 
the SDSS Data Release 7 (DR7; Abazajian et al. 2009). 
SDSS DR7 contains spectra of more than 300,000 stars. 
About 240,000 of them are from the Sloan Extension for 
Galactic Understanding and Exploration (SEGUE; Yanny 
et al. 2009), focused on the Galactic science. We select tar- 
get stars as those with line-of-sight extinctions Eb-v > 
0.3 and of spectra of signal-to-noise ratios S/N > 40. We 
restrict targets to objects with valid determinations of ef- 
fective temperature Toff, surface gravity log g and metal- 
licity [Fe/H] as given by the SEGUE Stellar Parameter 
Pipeline (SSPP, Lee et al. 2008a,b; AUende Prieto et al. 
2008; Lee et al. 2011; SmoUnski et al. 2011). Typical un- 
certainties of those parameters are 180 K, 0.25 dex and 
0.23 dex for T^s, log g and [Fe/H], respectively, dominated 
by systematic errors (Schlesinger et al. 2010; Smolinski et 
al. 2011). On the other hand, in the current work we are 
mainly interested in the relative ranking of stars in the T^g, 
log g and [Fe/H] parameter spaces, such that systematic 
uncertainties are not important. For robust template sub- 
traction, we further require the errors of T^e, log g and 
[Fe/H] to be smaller than 100 K, 0.2 dex and 0.1 dex, re- 
spectively. Values of Eb-v used here are from the extinction 
map of Schlegel et al. (1998). In total, 2,164 target stars are 
selected. Their spatial and stellar parameter distributions 
are displayed in Fig. 1. The targets come from spectroscopic 
plates that spread over a large sky area (—42° < 6 < -1-22°, 
8° ^ I ^ 210°), and involve stars of a wide range of spectral 
type, luminosity class and metallicity. Given the require- 
ments on S/N ratio and line-of-sight extinction, most tar- 
gets are from bright SEGUE plates towards low Galactic 
latitudes. Essentially, all the target stars are brighter than 
17 in g band. The requirements suggest that the majority 
of the target stars are from the Galactic disk, with the re- 
maining small fraction from the halo. This is consistent with 
the fact that most targets have [Fe/H] higher than —0.7, as 
shown in the bottom left panel of Fig. 1. 

For the control sample, we select stars of S/N ratios > 
20 and line-of-sight extinctions Eb-v < 0.03. The cuts on 
errors of Tefj , log g and [Fe/H] are the same as target stars. 
In total, 62,266 stars are selected. Given the large number 
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Figure 1. Properties of the sample of target stars. Top left: Spa- 
tial distribution in Galactic longitude I and latitude b. Top right: 
Tgff versus log g. Bottom left: T^g versus [Fe/H]. Bottom right: 
g — r versus g. 

of control stars, only those with the highest S/N ratios are 
used as templates, as described below. 

2.2 Method 

The idea of template subtraction method is straightforward. 

For each star in the target sample, stars of similar spec- 
tral types, luminosity classes and metallicities are selected 
from the control sample. Quantitatively, we require the dif- 
ferences in Teff, log g and [Fe/H] between a control and a 
target star to be smaller than 100 K, 0.5 dex and 0.2 dex, re- 
spectively. Secondly, after correcting for differences in radial 
velocities of control stars relative to the target star, spectra 
of the first one hundred control stars of highest S/N ratios 
are co-added to construct a template spectrum of very high 
S/N ratio. In cases where the number of control stars is less 
than one hundred, all control stars satisfying the above cri- 
teria are used. For the two DIB features of interest here, 
A5780 and A6283, a segment of the template spectrum of 
a width of 200 A centered around the feature of concern is 
scaled to the target spectrum, subtracted and then rectified 
using the continuum level determined from the scaled tem- 
plate spectrum. This certificated residual spectrum is devoid 
of stellar absorption features and enhances DIB features in 
cases of a good template matching. To account for the dif- 
ferent amount of extinction and other types of mismatch 
of the template and target spectra, the scaling factors of 
the segment of the template spectrum to target spectrum 
are fitted with a linear polynomial of constant slope after 
masking the DIB feature. Note that the strength of a DIB 
feature is very sensitive to the scaling factor at the wave- 
length of the feature, which is accurately estimated in the 
above fitting in cases of a good template matching. Given 
the relatively narrow wavelength span under consideration 
here, a continuum of constant slope is also used in fitting 
the continuum level of the scaled template spectrum. This 
continuum level is assumed to be equally applicable to the 
target spectrum. It is much easier to accurately determine 



the continuum level of the scaled template spectrum than 
that of the target, given the former generally has a very high 
S/N ratio that makes it easy to mask and exclude regions 
of stellar absorption features in the spectrum. Finally, we 
search for the DIB absorption feature around the expected 
wavelength in the rectified residual spectrum and measure 
its equivalent width using a line fitting method. 

The DIB A5780 contains two blended components, a 
narrow and a much broader one. The narrow and broad 
components have widths of 2.07 and 15.5 A centered at 
5780.59 and 5779.48 A, respectively (Jenniskens & Desert 
1994). Slightly different values, a width of 2.11 A centered at 
5780.48 A for the narrow component, were given by Hobbs ct 
al. (2008). Owing in part to the difficulty in reliably defining 
the contirmum level for the broad component, in the current 
work only the strengths of the narrow component will be 
considered, the DIB A5780 refers only to its narrow compo- 
nent. Given that the typical Doppler spread measured from 
interstellar atomic lines toward reddened stars {Eb-v < 0.8) 
is about 7.5 km s~^(Welty & Hobbs 2001) and the low spec- 
tral resolution of SDSS data, Doppler broadening is likely 
to be negligible. As a consequence, using a Gaussian profile 
of a fixed FWHM 3.85 A to fit and measure the strengths of 
the DIB A5780 is a good approximation. Here 3.85 A is the 
expected FWHM of the DIB A5780 after taking into account 
its intrinsic width and the SDSS instrumental broadening. 

The DIB A6283 also shows a relatively deep, narrow fea- 
ture superimposed on a rrmcli shallower, broader component. 
Examinations of the DIB A6283 observed with high resolving 
powers show that it exhibits a profile closely resembling a 
Lorentzian one, with its broad wings possibly contaminated 
by several weak overlapping DIB features (Jenniskens & 
Desert 1993; Drosback 2006). Hobbs ct al. (2008) estimated 
a FWHM of 4.77 A and a central wavelength of 6283.83 A for 
the DIB A6283 from the echelle spectra of the double-lined 
spectroscopic binary HD204820. Drosback (2006) obtained 
a FWHM of 4.0 ± 0.3 A from a larger sample. In this work, 
a fixed- width Voigt profile resulting from a convolution of a 
Lorentz profile of a FWHM 4.0 A and a Gaussian profile of 
a FWHM 3.49 A (the instrumental width of SDSS spectra), 
is used to fit and measure the strengths of DIB A6283. Note 
that the DIB A6283 feature is badly affected by telluric ab- 
sorption features. However, all released SDSS spectra have 
been corrected for telluric absorptions by dividing spectra 
of hot standard stars, so the measurements of DIB A6283 
in this work should not be much affected by telluric absorp- 
tions. 

Sample spectra and the resultant fits are illustrated in 
Fig. 2. Note in this Figure the spectra are plotted in vac- 
uum wavelengths. The standard formula given by Morton 
(1991) have been used to convert between air and vacuum 
wavelengths. 



3 RESULTS 

The DIBs AA5780, 6283 are detected in spectra of almost all 
target stars of a wide variety of types, from hot A to cool K 
types, dwarfs and giants, metal-rich and -poor stars. Dozens 
of target stars which lack obvious DIB features are nearby 
cool dwarfs and thus suffer little extinction [The extinction 
map of Schlegel et al. (1998) used to select target stars refers 
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Figure 2. The DIB features AA5780, 6283 detected in the SDSS spectra of six representative stars. Basic information of individual stars 
is hsted at the bottom of each row, including, from left to right, the SDSS spectroscopic target ID, spectral type, effective temperature, 
surface gravity, metallicity, line-of-sight extinction. Galactic longitude and latitude. For each row, segments of the target (black) and the 
scaled template (red) spectra centered at 5780 A and 6300 A are shown in the left two panels, respectively. The lines in red show fitted 
continuum levels of the scaled template spectra, assumed to be equally applicable to the target spectra. Those continua were used to 
rectified the target spectra after subtracting the scaled template ones. The residuals are plotted in the right two panels in black, with the 
Gaussian fitting to the DIB A5780 and Voigt fitting to the DIB A6283 over- plotted in red (see text for detail). Note that the DIB A5797 
is well detected in the spectra of the 1st, 3rd and 6th rows, and the DIB A6203 is well detected in the spectra of the 1st, 2nd and 4th 
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to the total extinction for a given sight-line]. Fig. 2 displays 
DIBs AA5780, 6283 detected in SDSS spectra of six represen- 
tative target stars. Basic information of individual example 
stars is listed in the Figure, including, from left to right, the 
SDSS spectroscopic target ID, spectral type, effective tem- 
perature, surface gravity, metallicity, line-of-sight extinction. 
Galactic longitude and latitude. The stellar parameters are 
those given by the SSPP pipeline, and extinction values are 
from Schlegel et al. (1998). As shown in the spectrum of a K7 
dwarf plotted in the 5th row, DIB A5780 suffers serious con- 
tamination from stellar absorption lines in cool and metal- 
rich stars, leading to some uncertainties in the identification 
and measurement of the feature. To verify the identification 
of the DIB A5780 feature in the spectrum of this K7 dwarf, 
we visually inspected individual spectra of its control stars 
and found that they agree with each other, confirming that 
the identification is real. Note that DIBs AA5797,6203 have 
been detected in more than half of the spectra. Several other 
DIB features, such as the DIBs AA5705, 5849, 6195, 6269, 
6308, 6317, 6376, 6379 possibly have also been detected in 
some of the spectra. 

For the DIBs AA5780, 6283, we have measured their 
depths and central wavelengths by Gaussian and Voigt pro- 
file fitting, respectively. The equivalent widths (EW) and 
radial velocities are also calculated. To derive the radial ve- 
locities, the rest-frame wavelengths of DIBs AA5780, 6283 
given by Hobbs et al. 2008 are used. Examples of the fit- 
ting are over-plotted in Fig. 2. In most cases, the target and 
template spectra match well and the fitting performs rea- 
sonably. In a few cases where there are not enough control 
stars to construct a matching template of high fidelity, the 
DIB features are weak, the spectra are of low S /N ratios or 
there are some defects in the spectra, the fitting procedure 
may deliver poor or even unreliable results. Nevertheless, we 
can still confirm the presence of the DIB features in most of 
those spectra by visual inspection. 

EW (5780) and £'W(6283) are very sensitive to the 
continuum levels, particularly for i?W (6283) because of its 
broad width. To estimate the uncertainties caused by con- 
tinuum determinations, we have re-fitted the DIBs AA5780, 
6283 with their continuum levels fixed to zero. For the 
DIB A5780, we find that values of EW{5780) measured 
with varying continuum levels are strongly correlated with 
(within an uncertainty of 0.03 A), but 10 per cent system- 
atically smaller than those measured with fixed ones. This 
is because the broad component of the DIB A5780 will gen- 
erally make the continuum level of its narrow component 
below zero when fitted with varying continuum levels. For 
the DIB A6283, values of £;W(6283) measured with varying 
continuum levels are also strong correlated with (within an 
uncertainty of 0.09 A), but 18 per cent systematically larger 
than those measured with fixed ones. This is because there 
are plenty of DIBs around the DIB A6283 within 100 A, 
including the relatively strong DIBs AA6195, 6203, 6269, 
6308, 6317, 6376, 6379. They wiU generally cause a slightly 
under-estimated continuum level of a target spectrum. The 
DIB A6283 feature is badly affected by telluric absorption 
features. All released SDSS spectra have been corrected for 
telluric absorptions. To check the uncertainties of SDSS tel- 
luric corrections and their effects on the EW (6283) measure- 
ments, we have selected a number of SDSS high S/N ratio 
spectra of low-extinction hot white dwarfs taken under dif- 
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Figure 3. Top left: Comparison of velocities measured from two 
DIB features A5780 and A6283, respectively. Top right: Com- 
parison of stellar and DIB A5780 velocities. Bottom left: Stellar 
(black) and DIB A5780 (red) velocities as a function of Galac- 
tic longitude. Bottom right: Stellar (black) and DIB A5780 (red) 
velocities as a function of Galactic latitude. To avoid crowding, 
only one in four randomly selected data points are plotted in the 
bottom left panel. 



ferent airmasses (1.0 - 1.2). The spectra have S/N ratios 
greater than 40 and Eb-v extinctions less than 0.03. Ex- 
aminations of those spectra show that the uncertainties of 
EW(6283) caused by any residuals in the telluric absorption 
corrections are likely to be well below 0.1 A. 

Detections of DIBs by Gaussian/ Voigt profile fitting 
may lead to false detections when the spectra are very noisy, 
or when stellar features are improperly subtracted. To avoid 
false detections, we have excluded a total number of 304 tar- 
gets that have either SVK(5780) < 0.082 A or EW {&28'S) < 
0.2 A in the following analysis. Although such a cut can not 
get rid of all false detections, the remaining ones, if any, are 
unlikely to affect the main results of the current work. 

The top left panel of Fig. 3 compares radial velocities 
deduced from the DIBs A5780 and A6283, respectively. It 
can be seen that the velocities deduced from DIB A6283 
are 16.0 km s~^systematically larger than those deduced 
from DIB A5780, which is likely caused by uncertainties of 
the adopted DIB rest-frame wavelengths from Hobbs et al. 
(2008). After accounting for the systematic differences, the 
two features give consistent results within an uncertainty of 
32 km s~^. Given that the DIB A5780 feature is narrower and 
more symmetric than the A6283 feature, the wavelengths of 
the DIB A5780 in the SDSS spectra are more easily to de- 
termine as well as its rest wavelength, so velocities deduced 
from the DIB A5780 are used to represent the DIB veloci- 
ties. The top right panel compares velocities deduced from 
the DIBs A5780 with velocities of the target stars from the 
SSPP pipeline. A weak correlation between the DIB and 
stellar velocities is seen. The correlation is due to the fact 
that most targets are disk stars and thus have kinematic 
characteristics similar to the absorbing interstellar clouds 
from which the DIBs originate. Also note compared to the 
wide range of stellar velocities, the span of the DIB velocities 
is much limited. This is expected considering the absorbing 
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Figure 4. Top left: EW{5780) versus Eb~v- Top right: _BVK(6283) versus Eb-v- Bottom panels: EW{57m) / EW{&2?,2,) versus 
-514^(5780). The black dots denote the targets of the current work. The red crosses donote the targets of the current work of spectral S/N 
ratios higher than 60. The red lines in the top panels represent regressions of the current work (using targets of spectral S/N ratios higher 
than 60). For comparison, the data and regressions of Friedman et al. (2011) are over-plotted in blue. In the bottom right panel, the red 
crosses are divided into 18 bins with a width of 0.03 A, the red squares in the red line indicate the median value for each bin. Note that 
the definitions of Eiy(5780) and SW(6283) in Friedman et al. (2011) and this work were very close, but EVF(5780) and EW(6283) in 
Friedman et al. were computed with direct integration. Also note that Eb-v values plotted in the top panels were all derived from the 
photometric colors of the target stars. 



clouds are confined to the Galactic thin disk and have more 
circular orbits than stars. That the DIB velocities follow the 
stellar values but with smaller scatters at a given longitude 
or latitude is clearly seen in the lower two panels of Fig. 3. 
Again, the trends are consistent with the fact that DIBs 
arise from the interstellar clouds in the Galactic thin disk. 

It is well known that to some extent strengths of in- 
dividual DIBs are correlated with each other, as well as 
the amount of intervening dust extinction. Using mea- 
surements from over one hundred stellar spectra of high 
S/N ratios and spectral resolution [R ~ 38,000), Fried- 
man et al. (2011) find that SW(5780) = 0.505 x Eb-v, 
SVF(6283) = 1.045 x Eb-v+ 0.080 and £^(5780) = 
0.431 X £1^(6283) - 0.012, with a correlation coefficient 
of 0.82, 0.82 and 0.96, respectively. Note that the defi- 
nitions of _EW(5780) and £;VK(6283) in Friedman et al. 
(2011) and this work were very close, but i?Vy(5780) and 
i?H^(6283) in Friedman et al. were computed with direct 
integration. Note also values of Eb-v of Friedman et al. 
(2011) were estimated from the observed and intrinsic col- 
ors. To investigate relations between £1^(5780), £I¥(6283) 



and Eb~v using measurements of stars in the target sam- 
ple of the current work, the largest ever assembled hith- 
erto, we plot _B'Py(5780) and EW{&29.2,) as a function of 
Eb-v and BW" (5870) /fiW (6283) against £1^(5780) in 
Fig. 4. For comparison, data and relations from Friedman 
et al. (2011) are over-plotted. 

Given that values of Eb-v from the extinction map of 
Schlegel et al. (1998) suffer from the following limitations: 

a) The map gives the total amount of extinction integrated 
along a given line-of-sight to infinite, and consequently the 
value is an upper limit of the real one for a local disk star; 

b) The map fails at low Galactic latitudes (|b| ^ 5°); and 

c) The map has a limited spatial resolution about 6 ar- 
cmin, whereas the extinction may well vary below such a 
scale, we did not use the values of Eb-v from Schlegel et 
al. in Fig. 4. The Eb-v value of a target used here was de- 
rived directly from the difference of its observed and intrinsic 
g — r colors. The intrinsic g — r color was chosen to be the 
median value of the dereddened g — r colors of its control 
stars. Here Ag/EB-v = 3.303 and Ar/EB-v = 2.285 at a 
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total-to-selective extinction ratio of 3.1 were used (Schlafly 
& Finkbeiner 2011). 

The variations of £^(5780) and EW{6283) as a func- 
tion of Eb-v deduced from the current work are shown in 
the top left and right panels of Fig. 4, respectively. The black 
dots denote the targets of the current work. The red crosses 
denote the targets of the current work of spectral S/N ra- 
tios higher than 60. Linear regressions passing through the 
origin were performed on the red crosses in the top panels, 
with each data point carrying equal weight. The relations 
obtained, i.e., EW{b7%0) = 0.61 x Eb-v and £^(6283) 
= 1.26 X Eb-v, are plotted in red lines. The blue crosses 
and lines denote the data of and the relations obtained by 
Friedman et al. (2011), respectively. The slopes of regres- 
sions of this work are about 20 per cent larger than those of 
Friedman et al. (2011). However, considering the large scat- 
ters in both data, this work yields consistent results with 
Friedman et al. (2011). 

Values of the £;W(5780)/SW^(6283) ratios as deduced 
for the current target sample stars are plotted against 
iJW(5780) in the bottom panels of Fig. 4. The symbols are 
the same as in the top panels. In the bottom right panel, red 
crosses are divided into 18 bins with a width of 0.03 A, red 
squares in the red line indicate the median value for each bin. 
The data in this work yield a median value of 0.48, consis- 
tent with the result of ~ 0.42 of Friedman et al. (2011). The 
data show large scatters, possibly due to the limited S/N 
ratios of the target spectra. The scatters of data deduced 
from spectra of S/N ratios higher than 60 are comparable 
to those of Friedman et al. (2011), implying measurements of 
the current work are robust. The red line in the bottom right 
panel seems to indicate that (5780) /BW (6283) first in- 
creases with EW{h79>Q) and then saturates for £1^(5780) 
above some value between 0.2 and 0.3 A. This interesting 
trend is supported by the data of Friedman et al. (2011). 
Friedman et al. attribute the behaviors of the strength of 
the DIB A5780 relative to that of A6283 to measurement er- 
rors or a threshold effect, where one DIB cannot form until a 
significant amount of another DIB is present. Additionally, 
it could be interpreted that the fractions of carrier of indi- 
vidual DIBs depend on certain local environment factor(s) 
in an equilibrium state. For example, Vos et al. (2011) find 
that the EW {^ISQ) / EW ratio depends on cloud den- 
sity as well as exposure to the interstellar radiation field. 

Taking into account the uncertainties of _EH^(5780), 
i?W(6283) and Eb-v, the results plotted in Fig. 4 are con- 
sistent with the data obtained by Friedman et al. (2011), 
particularly for those deduced from spectra of S/N ratios 
higher than 60. Both the data of the current work and Fried- 
man et al. (2011) show large scatters, suggesting that a sig- 
nificant fraction of the scatters are intrinsic. We will investi- 
gate whether the properties of DIBs vary with Galactic posi- 
tion and environment in a separate paper, using Eb-v val- 
ues and distances of stars determined directly from photo- 
metric and spectroscopic observations. 



4 DIBS IN THE LAMOST COMMISSIONING 
SPECTRA 

The Guoshoujing Telescope (the Large Sky Area Multi- 
Object Fiber Spectroscopic Telescope LAMOST; Wang 



et al. 1996; Su et al. 1998; Xing et al. 1998; Zhao 
2000; Cui et al. 2004; Zhu ej al 2 006; Cui et al. 2012; 
c.f. http:/7www.lamost.org/website/en/ ), is an innovative 
quasi-meridian reflecting Schmidt Telescope that is capable 
of recording spectra of up to 4,000 objects simultaneously 
in a field of view of 5° in diameter. Following a two-year 
commissioning initiated in September, 2009 and a one-year 
pilot survey initiated in October, 2011, the LAMOST for- 
mal Galactic and extragalactic surveys are expected to com- 
mence in fall of 2012 (Zhao et al. 2012). 

Given the similar spectral resolving powers of SDSS and 
LAMOST, DIBs should be also detectable in the LAMOST 
spectra. Given that the data reduction pipelines of LAM- 
OST are still under development and reliable stellar pa- 
rameter estimates are not available yet, the template sub- 
traction method currently does not work for the LAMOST 
data. However, taking advantage of the smooth stellar con- 
tinua of hot early-type stars, DIBs are also easily identi- 
fied in the LAMOST low-resolution spectra. As an exam- 
ple, at least nine DIB features (AA4429, 5780, 5797, 6203, 
6283, 6532, 6613, 6993, 8621) have been detected in the spec- 
tra of a LAMOST target (Fig. 5). Note the prominence of 
DIBs AA5780, 6283. The target was observed on December 
21, 2009 with the half-slit (1.65 arcsec) mode (R ~ 2,000) 
during the commissioning phase of LAMOST. Three expo- 
sures, each of 900s, were obtained. The data were reduced us- 
ing the LAMOST 2D pipeline (Luo, Zhang & Zhao 2004) fol- 
lowing the standard procedure, including bias subtraction, 
cosmic-ray removal, ID spectral extraction, flat-fielding and 
wavelength calibration. Sky subtraction was not performed. 
This did not affect the identification of DIB features given 
the weakness of sky continuum level compared to that of the 
target. Flux calibration was also not necessary. It is worth 
mentioning that the target is an emission line star first dis- 
covered by Kohoutek & Wehmeyer (1999). Its broad and 
prominent Hq emission line is clearly seen in the bottom 
left panel of Fig. 5. 



5 DISCUSSION 

Using a template subtraction method, we have successfully 
identified the DIBs AA5780, 6283 in the SDSS low-resolution 
spectra of a sample of about 2,000 stars and measured their 
equivalent widths and velocities. The sample is by far the 
largest of stars with DIBs detected and measured. The sam- 
ple stars are distributed over a large sky area and cover a 
wide range of stellar parameters of effective temperature, 
surface gravity and metallicity, confirming that the carriers 
of DIBs are ubiquitous in the diffuse ISM. We use the sample 
to investigate the relations between the equivalent widths of 
the DIBs AA5780, 6283 and sight-line extinctions, and find 
results consistent with previous works. The data however 
show large scatters. Some of the scatters are due to the mea- 
surement uncertainties, especially for spectra of limited S/N 
ratios. A significant fraction of the scatters are probably in- 
trinsic - the data of Friedman et al. (2011) also show signif- 
icant scattering. Identifying potential environmental factors 
that may lead to the complexities of observed properties of 
DIBs is pivotal to reveal the nature of DIBs. To achieve 
such a goal, accurate measurements of DIBs towards a large 
number of sight-lines probing a variety of environments with 
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Figure 5. Nine DIB features have been detected in the LAMOST commissioning spectra of an emission line star. The spectra have 
been normaUzed. The different colors in each panel refer to individual spectra from three separate exposures. 



low-cost spectra of low and intermediate resolution as illus- 
trated in the current work should be an effective approach. 

Apart from the limited S/N ratios of spectra, other 
sources of uncertainties in measuring the DIBs AA5780, 6283 
in this work mainly include a mismatch between a target 
and template spectrum, continuum determinations and pro- 
file fitting of the DIB features with a single Gaussian/Voigt. 
A mismatch of a target and template spectrum may occur 
if the target is of a rare spectral type, the stellar parame- 
ters of the target star are in significant error or there are 
some defects (e.g., residual telluric absorption features) in 
its spectrum. Such cases can in principle be flagged out and 
excluded from the sample in the future. A mismatch can also 
happen if the distribution of stellar parameters of the control 
sample is biased. Such cases can be avoided by more care- 
ful selections of control stars. To reduce the uncertainties of 
continuum levels, sophisticated ways to define stellar con- 
tinuum are needed. To reduce the uncertainties introduced 
by a single Gaussian/Voigt profile fitting, one can explore 
deblending the DIB A5780 with two Gaussians and fitting 
the DIB A6283 with measured profiles. It might be possi- 
ble to carry out the fitting with the central wavelengths of 
individual components of the DIBs tied together, reducing 
the number of degrees of freedom of the fitting function. For 



spectra of high S/N ratios, it might be feasible to define 
indexes of DIBs and measure them by direct integration. 

It is worth mentioning that hot stars have few strong 
absorption lines blending with DIBs. In such cases, template 
spectra are not even needed to identify DIBs. As illustrated 
by the example presented in the previous Section, in addi- 
tion to the strong DIBs AA5780, 6283, seven other DIBs can 
be easily identified in low-resolution spectra of such a hot 
star. The strengths of the DIB features can also be measured 
by direct integration given sufficient high S/N ratios. Given 
that the light of a star-burst galaxy is dominated by con- 
tribution from luminous hot stars, searching for the strong 
DIBs (e.g., DIBs AA5780, 6283) in their SDSS spectra might 
be feasible. 

Searching and measuring DIBs in low-resolution spec- 
tra, now widely available for a huge number of objects, have 
obviously the potential to tackle a number of problems of 
the ISM. For example, they can be used to probe the distri- 
bution and properties of the DIBs, estimate the extinction 
towards individual stars, trace the amount of ISM absorp- 
tion and search for extremely metal-poor stars. The data 
will help unveil the nature of DIBs. Large scale spectro- 
scopic surveys, such as the on-going RAVE and SDSS, and 
the up-coming LAMOST Galactic spectral surveys will yield 
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detections of DIBs (e.g., AA4429, 5780, 5797, 6203, 6283, 
6613, 6993, 8621) in hundreds of thousands stars of a va- 
riety of spectral types, located in different parts and envi- 
ronments of the Milky Way. The huge database will provide 
an unprecedented opportunity to study the demographical 
distribution and nature of DIBs as well as using DIBs to 
probe the distribution and properties of the ISM and the 
dust extinction. Given the large sample, targets exhibiting 
particularly interesting features can be selected and followed 
with high-resolution spectroscopic observations. 

Strengths of DIBs measured with low-resolution spec- 
troscopy can be used to estimate extinction towards a large 
number of individual stars targeted in large scale spec- 
troscopic surveys. Extinction is one of the most funda- 
mental stellar parameters and affects the determinations 
of others (i.e. effective temperature, metallicity and dis- 
tance). As mentioned above, extinction values given by 
Schlegel et al. (1998) have several limitations, especially for 
the low Galactic latitude regions. From the relations be- 
tween the strengths of DIBs and the line-of-sight color ex- 
cess, strengths of DIBs measured with low-resolution spec- 
troscopy provide an independent way to estimate extinc- 
tion of individual stars. This is especially useful for the up- 
coming LAMOST Digital Sky Survey of the Galactic Anti- 
center (DSS-GAC; Liu et al. 2012, in prep). The DSS-GAC 
project aims to obtain spectra for a statistically complete 
sample of over three million stars down to 18.5 limiting 
magnitude in r- band, distributed in a contiguous area of 
3438 square degree (-30° < 6 < 4-30°, 150° < / < 210°) 
and sampling a significant volume of the Galactic thin/thick 
disks, halo and their interfaces. The project will target ap- 
proximately 1,000 stars per square degree outside the Galac- 
tic plane (|6| 3.5°) and twice that within the Galactic 
plane. Although extinction towards Galactic anti-center is 
much lower than towards the center, there are regions of 
high extinction where the extinction map of Schlegel et al. 
(1998) is inapplicable. On the other hand, information of 
extinction is required to select spectral standard stars and 
perform flux calibration. Information of extinction is also 
required to interpret the target selection algorithms of the 
survey and derive basic stellar parameters and the selection 
function of the survey. The possibility of estimating the ex- 
tinction to individual star via the observed DIB strengths, 
unaffected by flux calibration, should be very helpful to the 
DSS-GAC project. In the meantime, the DSS-GAC project 
will deliver a huge database to study DIBs. The database 
will generate a three-dimensional extinction map towards 
the Galactic anti-center, with distances of individual stars 
estimated from either a spectro-photometric analysis or di- 
rectly from the parallax measurements to be delivered by the 
forth-coming GAIA mission (Ferryman et al. 2001). Finally, 
given the tight correlations between DIB strengths (espe- 
cially the DIB A5780; Friedman et al. 2011) and the neutral 
hydrogen column density, the database can be further used 
to derive the three-dimensional distribution of the neutral 
hydrogen, probe the properties, structure and kinematics of 
the ISM, and constrain the peculiar motion of the sun with 
respect to the Local Standard of Rest. 

It is possible that detections of DIBs with low-resolution 
spectroscopy may provide a way to correct for contamination 
of the ISM Ca II absorption when searching for extremely 
metal-poor stars with low-resolution spectroscopic surveys. 



Discovery of extremely metal-poor stars is of great impor- 
tance in understanding the nature and the initial mass func- 
tion of the flrst stars, the processes of element production 
during the early evolution phase of the Milky Way and the 
metallicity distribution function of the Galactic halo (Beers 
& Christlieb, 2005). Low- resolution spectroscopic surveys 
have been playing an important role in searching for very 
metal-poor stars (Beers & Caxollo 2008; Beers et al 2009; 
Fulbright et al. 2010). Under low spectral resolution, the 
stellar Ca II absorption lines blend with absorption aris- 
ing from intervening ISM, leading to overestimated stellar 
metallicities. This problem becomes more severe as stars get 
more metal-poor. Strengths of certain DIBs might provide 
a way to correct for the ISM Ca II absorption. Using high- 
resolution spectra of a sample of 35 stars, Galazutdinov et 
al. (2004) show that the strengths of some DIBs correlate 
well with the interstellar K I lines but to a much lower de- 
gree with the Ca II lines. From a sample of several hundred 
low-resolution spectra of 0-type stars, Penades Ordaz et al. 
(2011) flnd a moderately low correlation coefflcient between 
the Ca II K line absorption and the extinction with a pecu- 
liar spatial distribution that they ascribe to a relationship 
between absorption saturation and velocity profile of the 
Ca II K line. Nevertheless, this is a problem worth further 
investigations. 

Finally, we would like to point out that the template 
subtraction technique used in the current work can be 
adapted to tackle other problems, such as searching for pecu- 
liar stars of unusual characteristics in a large spectroscopic 
database where templates of "normal" stars can be easily 
constructed. Examples of such peculiar objects include binsr 
ries, emission line stars and stars of peculiar chemical com- 
position. 
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